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Abstract 

We introduce a new determinantal method to count cycle systems in a directed graph that gen¬ 
eralizes Gessel and Viennot’s determinantal method on path systems. The method gives new insight 
into the enumeration of domino tilings of Aztec diamonds, Aztec pillows, and related regions. 


1 Introduction 

In this article, we present an analogue of the Gessel-Viennot method for counting cycle systems on a type 
of directed graph we call a hamburger graph. A hamburger graph H is made up of two acyclic graphs Gi 
and G2 and a connecting edge set £3 with the following properties. The graph Gi has k distinguished 
vertices {ui,... ,Vk} with directed paths from Vi to Vj only if i < j. The graph G2 has k distinguished 
vertices {wk+i, ■ ■ ■ ,W2k} with directed paths from Wi to Wj only if i > j. The edge set £3 connects 
each vertex Vi to vertex Wk+i and vice versa. (See Figure ^ for a visualization.) Hamburger graphs arise 
naturally in the study of Aztec diamonds, as explained in Section 0 

The Gessel Viennot method is a determinantal method to count path systems in an acyclic directed 
graph G with k sources S\,... ,Sk and fc sinks ..., tk- A path system £ is a collection of k vertex-disjoint 
paths, each one directed from s, ; to for some permutation a £ Sk (where Sk is the symmetric group 
on k elements). Call a path system V positive if the sign of this permutation a satisfies sgn(cr) = +1 
and negative if sgn(cr) = —1. Let p + be the number of positive path systems and p~ be the number of 
negative path systems. 


j ] j ] j ] 


Figure 1 : A hamburger graph 
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Figure 2: A simple hamburger graph H 


Corresponding to this graph G is a k x k matrix A = (ciij), where is the number of paths from s, 
to tj in G. The result of Gessel and Viennot states that det A = p + — p~. The Gessel-Viennot method 
was introduced in J3i i] , and has its roots in works by Karlin and McGregor [5| and Lindstrom fSJ ■ A 
nice exposition of the method and applications is given in the article by Aigner jT| . 

This article concerns a similar determinantal method for counting cycle systems in a hamburger graph 
H. A cycle system C is a collection of vertex-disjoint directed cycles in H. Let l be the number of edges 
in C that travel from Gi to G\ and let m be the number of cycles in C. Call a cycle system positive if 
(— l) l + m = +1 and negative if (—l) ,+m = —1. Let c + be the number of positive cycle systems and c~ 
be the number of negative cycle systems. Corresponding to each hamburger graph H is a 2k x 2k block 
matrix Mh of the form 


M h = 


A I k 
-h B 


where in the upper triangular matrix A = (aij), a%j is the number of paths from Vi to v 3 in Gi and in 
the lower triangular matrix B = (bij), bij is the number of paths from Wk+i to Wk+j in G 2 . This matrix 
Mh is referred to as a hamburger matrix. 


Theorem 1.1 (The Hamburger Theorem). If H is a hamburger graph , then det Mh = c + — c . 

A hamburger graph H is called strongly planar if there is a planar embedding of H that sends Vj to 
(■ i , 1) and w^+i to ( i , —1) for all 1 < i < k, and keeps edges of E\ in the half-space y > 1 and edges of E 2 
in the half-space y < —1. This definition suggests that Gi and G 2 are “relatively” planar in H, a stronger 
condition than planarity of H. Notice that when H is strongly planar, each cycle must use exactly one 
edge from G 2 to Gi. Hence, the sign of every cycle system is +1. This implies the following corollary. 


Corollary 1.2. If H is a strongly planar hamburger graph, det Mh = c + . 


The following simple example serves to guide us. Consider the two graphs Gi = (Vi,E\) and G 2 = 
(V 2 ,E 2 ), where V\ = {vi,v 2 ,V 3 }, V 2 = {wi,W 5 ,we}, E\ = (ni —> v 2 ,v 2 —*• v%,v\ — > V 3 }, and E 2 = {we —** 
W 5 ,we W 4 ,We —> W 4 }. Our hamburger graph H will be the union of Gi, G 2 , and the edge set E 3 
consisting of edges e* : Vi —> Wk+i and e' : Wk+i —> Vi- In this example, k = 3 and H is strongly planar. 
Figure |3 gives a graphical representation of H. 

In this example, the hamburger matrix Mh equals 


M h 


1 1 2 10 0 
0 1 10 10 

0 0 1 0 0 1 
-10 0 10 0 
0 -1 0 110 

0 0-1211 


The determinant of Mh is 17, corresponding to the seventeen cycle systems (each with sign +1) in Figure 

El 
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Figure 3: The seventeen cycle systems for the hamburger graph in Figure |3 


The graph that inspired the definition of a hamburger graph comes from the work of Brualdi and 
Kirkland [2j , in which they give a new proof that the number of domino tilings of the Aztec diamond is 
2"( n+ i)/ 2 . An Aztec diamond , denoted by AD n , is the union of the 2 n(n + 1) unit squares with integral 
vertices (x,y) such that \x\ + \y\ < n + 1. An Aztec pillow, as it was initially presented in M, is also a 
rotationally symmetric region in the plane. On the top left boundary, however, the steps are composed 
of three squares to the right for every square up. Another definition is that Aztec pillows are the union of 
the unit squares with integral vertices (x,y) such that |a; + y| < n + 1 and |3y —x| < n + 3. As with Aztec 
diamonds, we denote the Aztec pillow with 2 n squares in each of the central rows by AP n . In Section 0 
we extend the notion of Aztec pillows having steps of length 3 to “odd pillows”—those that have steps 
that are of a constant odd length. The integral vertices (x, y ) of the unit squares in g-pillows for q odd 
satisfy \x + y\ < n + 1 and \qy — x\ < n + q. 

We introduce the idea of a generalized Aztec pillow, where the steps on all diagonals are of possibly 
different odd lengths. More specifically, a generalized Aztec pillow is a horizontally convex and vertically 
convex region such that the steps both up and down in each diagonal have an odd number of squares 
horizontally for every one square vertically. A key fact that we will use is that any generalized Aztec 
pillow can be recovered from a large enough Aztec diamond by the placement of horizontal dominoes. 
See Figure 01 for examples of an Aztec diamond, an Aztec pillow, and a generalized Aztec pillow. 

Brualdi and Kirkland prove the formula for the number of domino tilings of an Aztec diamond by 
creating an associated digraph and counting its cycle systems, manipulating the digraph’s associated 
Kasteleyn- Percus matrix of order n(n + 1). To learn about Kasteleyn theory and Kasteleyn-Percus 
matrices, start with Kasteleyn’s 1961 work 7] and Percus’s 1963 work 0. The Hamburger Theorem 
proves that we can count the number of domino tilings of an Aztec diamond with a much smaller 
determinant, of order 2 n. An analogous reduction in determinant size (from order 0(n 2 ) to order 0(n)) 
occurs for all regions to which this theorem applies, including generalized Aztec pillows. In addition, 
whereas Kasteleyn theory applies only to planar graphs, there is no planarity restriction for hamburger 
graphs. For this reason, the Hamburger Theorem gives a new counting method for cycle systems in some 
non-planar graphs. 


THE ELECTRONIC JOURNAL OF COMBINATORICS 12 (2005), #R00 


3 


Figure 4: Examples of an Aztec diamond, an Aztec pillow, and a generalized Aztec pillow 


In Section [2 we present an overview of the proof of the Hamburger Theorem, including the key 
lemmas involved. The necessary machinery is built up in Section 0 to complete the proof in Section 
01 Section 0 presents applications of the Hamburger Theorem to Aztec diamonds, Aztec pillows, and 
generalized Aztec pillows. Section 01 concludes with a counterexample to the most natural generalization 
of the Hamburger Theorem and an extension of Propp’s Conjecture on Aztec pillows. 

2 Outline of the Proof of the Hamburger Theorem 

2.1 The Hamburger Theorem 

Like the proof of the Gessel-Viennot method, the proof of the Hamburger Theorem hinges on cancel¬ 
lation of terms in the permutation expansion of the determinant of Mh- In the proof, we must allow 
closed directed walks in addition to cycles. We must also allow walk systems , arbitrary collections of 
closed directed walks, since they can and will appear in the permutation expansion of the hamburger 
determinant. We call a walk system simple if the set of walks visits no vertex more than once. We call a 
cycle of the form c : Vi —> Wi+k —> Vi a 2 -cycle. 

Each signed term in the permutation expansion of the hamburger determinant is the contribution 
of many signed walk systems W. Walk systems that are not cycle systems will all cancel out in the 
determinant expansion. We will show this in two steps. We start by considering walk systems that are 
not simple. If this is the case, one of the two following properties MAY hold. 

Property 1. The walk system contains a walk that has a self-intersection. 

Property 2. The walk system has two intersecting walks, neither of which is a 2-cycle. 

The following lemma shows that the contributions of walk systems satisfying either of these two 
properties cancel in the permutation expansion of the determinant of Mh- 

Lemma 2.1. The set of all walk systems W that satisfy either Property 1 or Property 2 can he parti¬ 
tioned into equivalence classes, each of which contributes a net zero to the permutation expansion of the 
determinant of Mh ■ 

The proof of Lemma l‘2.1 l oses a generalized involution principle. Walk systems cancel in families based 
on the their “first” intersection point. 

The remainder of the cancellation in the determinant expansion is based on the concept of a minimal 
walk system; we motivate this definition by asking the following questions. What kind of walk systems 
does the permutation expansion of the hamburger determinant generate, and how is this different from 
our original notion of cycle systems that we wanted to count in the introduction? The key difference is 
that the same collection of walks can be generated by multiple terms in the determinantal expansion of 
Mh ; whereas, we would only want to count it once as a cycle system. This redundancy arises when the 
walk visits three distinguished vertices in G\ without passing via G 2 or vice versa. We illustrate this 
notion with the following example. 
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Consider the second cycle system in the third row of Figure^ consisting of one solitary directed cycle. 
Since this cycle visits vertices Vi, v 2 , V3, iuq , and W4 in that order, it contributes a non-zero weight in 
the permutation expansion of the determinant corresponding to the term (12364) in Sq . Notice that this 
cycle also contributes a non-zero weight in the permutation expansion of the determinant corresponding 
to the term (1364). We see this since our cycle follows a path from v\ to V 3 (by way of v 2 ), returning to v\ 
via We and W4. We must deal with this ambiguity. We introduce the idea of a minimal permutation cycle, 
one which does not include more than two successive entries with values between 1 and k or between 
k + 1 and 2k. We see that (1364) is minimal while (12364) is not. 

We notice that walk systems arise from permutations, so it is natural to think of a walk system first 
as a permutation, and afterward a collection of walks that “follow” the permutation. This is the idea of 
a walk system-permutation pair (or WSP-pair for short) that is presented in Section EP1 From the idea 
of a minimal permutation cycle, we define a minimal walk to have as its base permutation a minimal 
permutation cycle, and a minimal walk system to be composed of only minimal walks. Since our original 
goal was to count “cycle systems” in a directed graph, we realize we need to be precise and instead count 
“simple minimal walk systems”. This leads to the second part of the proof of the Hamburger Theorem. 

Given a walk system that satisfies neither Property 1 nor Property 2 and is neither simple nor minimal, 
at least one of the two following properties MUST hold. 

Property 3. The walk system has two intersecting walks, one of which is a 2-cycle. 

Property 4. The walk system is not minimal. 

The following lemma shows that the contributions of walk systems satisfying either of these new 
properties cancel in the permutation expansion of the determinant of Mh- 

Lemma 2.2. The set of all walk systems W that satisfy neither Property 1 nor Property 2 and that satisfy 
either Property 3 or Property 4 can be partitioned into equivalence classes, each of which contributes a 
net zero to the permutation expansion of the determinant of Mh ■ 

The proof of Lemma E21 is also based on involutions. Walk systems cancel in families built from an 
index set containing the set of all 2-cycle intersections and non-minimalities. 

If a walk system satisfies none of the conditions of Properties 1 through 4, then it is indeed a simple 
minimal walk system, or in other words, a cycle system. The cancellation from the above sets of families 
gives that only cycle systems contribute to the permutation expansion of the determinant of Mh- This 
contribution is the signed weight of each cycle system, so the determinant of Mh exactly equals c + — c~. 
Theorem o follows from Lemmas EH and 12.21 in Section @1 • 


2.2 The Weighted Hamburger Theorem 

There is also a weighted version of the Hamburger Theorem, and it will be under this generalization 
that Lemmas EU and IQ are proved. We allow weights wt(e) on the edges of the hamburger graph; 
the simplest weighting, which counts the number of cycle systems, assigns wt(e) = 1. We require that 
wt(ej)wt(e() = 1 for all 2 < i < k — 1, but we do not require this condition for i = 1 nor for i = k. Define 
the 2 k x 2 k weighted hamburger matrix Mh to be the block matrix 


M h 


A D\ 
-D 2 B 


(1) 


In the upper-triangular k x k matrix A = (), a»j is the sum of the products of the weights of edges 
over all paths from Vi to Vj in Gi. In the lower triangular k x k matrix B = (6^), bij is the sum of the 
products of the weights of edges over all paths from Wk+i to Wk+j in G 2 . The diagonal k x k matrix D\ 
has as its entries da = wt(e*) and the diagonal k x k matrix D 2 has as its entries da = wt(e(). Note that 
when the weights of the edges in E$ are all 1, these matrices satisfy D\ = D 2 = /&. 
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We wish to count vertex-disjoint unions of weighted cycles in H. In any hamburger graph H , there 
are two possible types of cycle. There are k 2-cycles 

c : Vi u>k+i Vi 

and many more general cycles that alternate between G\ and G 2 . We can think of a general cycle as a 
path Pi in G\ connected by an edge eqi S £3 to a path Qi in G 2 , which in turn connects to a path P 2 
in Gi by an edge e( 2 , continuing in this fashion until arriving at a final path Qi in G 2 whose terminal 
vertex is adjacent to the initial vertex of P\. We write 


Pi 


ei.i 


Q 1 


P 2 


Pi 


Qi- 


For each cycle c, we define the weight wt(c) of c to be the product of the weights of all edges traversed 
by c: 

wt(c) = JJ wt(e). 

e€=c 

We define a weighted cycle system to be a collection C of m vertex-disjoint cycles. We again define the 
sign of a weighted cycle system to be sgn(C) = (—l) i+m , where l is the total number of edges from G 2 to 
Gi in C. We say that a weighted cycle system C is positive if sgn(C) = +1 and negative if sgn(C) = —1. 
For a hamburger graph H, let c + be the sum of the weights of positive weighted cycle systems, and let 
be the sum of the weights of negative weighted cycle systems. 

Theorem 2.3 (The weighted Hamburger Theorem). The determinant of the weighted hamburger 
matrix Mh equals c + — c~. 

As above, Theorem l2.3l follows from Lemmas l 2 .ll and P 2 . 2 l The proofs will be presented after developing 
the following necessary machinery. 


3 Additional Definitions 

3.1 Edge Cycles and Permutation Cycles 

In the proof of the Hamburger Theorem, there are two distinct mathematical objects that have the name 
“cycle”. We have already mentioned the type of cycle that appears in graph theory. There, a (simple) 
cycle in a directed graph is a closed directed path with no repeated vertices. 

Secondly, there is a notion of cycle when we talk about permutations. If a £ S n is a permutation, we 
can write a as the product of disjoint cycles a = X 1 X 2 '' ■ Xt■ 

To distinguish between these two types of cycles when confusion is possible, we call the former kind an 
edge cycle and the latter kind a permutation cycle. Notationally, we use Roman letters when discussing 
edge cycles and Greek letters when discussing permutation cycles. 

3.2 Permutation Expansion of the Determinant 

We recall that the permutation expansion of the determinant of an n x n matrix M = () is the 
expansion of the determinant as 


det M = ^ (sgn ct)to 1jCT(1) • • • . (2) 

crG Sn 

We will be considering non-zero terms in the permutation expansion of the determinant of the hamburger 
matrix Mh■ Because of the special block form of the hamburger matrix in Equation 0 , the permutations 
(7 that make non-zero contributions to this sum are products of disjoint cycles of either of two forms—the 
simple transposition 

X = (¥>11 wn) 
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or the general permutation cycle 

X=0Pll Pl2 ••• Ull U> 12 ■■■ U lvi ¥>21 . v\l ••• V\v x )- (3) 

In the first case, wn = ipu + k. In the second case, 1 < <p LK < k, k + 1 < u LK < 2k, (p LK < tp hK +i, and 
u> lK > <4 >i,k +i for all 1 < t < A and relevant k. The block matrix form also implies that (p L/J , L + k = to L i, 
ui IVl — k = and ui\ Ux — k = <pu. So that this permutation cycle is in standard form, we make sure 

that <^>n = min t)K <p LK . In order to refer to this value later, we define a function $ by $(%) = <pn- Each 
value 1 < ip t < k or k + 1 < to L < 2k appears at most once for any er £ S 2 k- 

We call a permutation cycle x minimal if it is a transposition or if /x t = v L = 2 for all i. This condition 
implies that we can write our general permutation cycles \ i n the form 

X = (v>ll ^12 Wn W12 ¥>21 • ■ • <*Al U\ 2 ), (4) 

with the same conditions as before. We call a permutation cr = Xi ■ • ■ XV minimal if each of its cycles Xt 
is minimal. 


3.3 Walks Associated to a Permutation 

To each permutation x G S 2 k, we can associate one or more walks c x in H. 

If x is the transposition x = (^n uu), then we associate the 2-cycle c x : —> w uil —> v vil to X- 

To any permutation x that is n °t a transposition, we can associate multiple walks c x by gluing together 
paths that follow x i n the following way. If x has the form of Equation ©, then for each 1 < * < A, 
let Pi be any path in G± that visits each of the vertices v Vil , v Vi2 , all the way through v Vi in order. 
Similarly, let Qi be any path in G 2 that visits each of the vertices w Uil , w Wi2) through w^. in order. For 
each choice of paths Pi and Qi, we have an additional possibility for the walk c x ; we can set 


c x : P\ 





Qx- 


(5) 


See Figure 0 for the choices of C(i 2 364 ) in the hamburger graph presented in Figure 0 We call A the 
number of P-paths in c x . The function <f>, defined in the previous section, defines a partial ordering 
on walks in a walk system—we say that the associated walk c x comes before the associated walk c x ' if 
<h(y) < $(x'). We call this the initial term order. As in Section 0H we define the weight of a walk c x to 
be the product of the weights of all edges traversed by c x . 


3.4 Walk System-Permutation Pairs 

We defined walk systems in Section 2, but we will see that the proof of Theorem o requires us to 
think of walk systems first as a permutation and second as a collection of walks determined by the 
permutation. We will see that for cycle systems as presented initially, signs and weights are not changed 
by this recharacterization. 

If H is a hamburger graph with k pairs of distinguished vertices, we define a walk system-permutation 
pair as follows. 

Definition 3.1. A walk system-permutation pair (or WSP-pair for short) is a pair (W, cr), where er £ S 2 k 
is a permutation and W is a collection of walks c £ W with the following property: if the disjoint cycle 
representation of a is er = ''' Xt , then W is a collection of r walks c X l , for 1 < u < r, such that c Xl is 

a walk associated to the permutation cycle y t . 

We define the weight of a WSP-pair (W, cr) to be the product of the weights of the associated walks 
c x £W. 

The definition of WSP-pair implies that each permutation er yields many collections of walks W, that 
collections of walks W may be associated to many permutations er, but that any simple walk system W 
corresponds to one and only one minimal permutation cr m . This is because, given any path as in Equation 
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Figure 5: A permutation cycle x and the two walks in H associated to x 


©, we can read off the initial and terminal vertices of each Pi and Qi in order, producing a well-defined 
permutation cycle a m . We define a WSP-pair (>V, cr) to be minimal if cr is a minimal permutation. 

For a WSP-pair (W, cr), where a = xi ' • ‘Xt, we define the sign of the WSP-pair, sgn(W, a), to be 
(—l) z sgn(cr), where A x is the number of P-paths in c x and where l = X)c x ew ^x- Alternatively, we could 
consider the sign of (W, cr) to be the product of the signs of its associated walks c x , where the sign of c x 
is sgn(c x ) = (—l) A *sgn(x). We say that a WSP-pair (W, cr) is positive if sgn(W, a) = +1 and is negative 
if sgn(>V, cr) = —1. 

Note that if (W,cr) is a minimal WSP-pair, then sgn(c x ) = +1 for a transposition y and sgn(c x ) = 
(—1) A+1 if x is °f the form in Equation 0. In particular, when (W, cr) is minimal and simple, its sign 
and weight is consistent with the definition given in the introduction. 


4 Proof of the Hamburger Theorem 

As mentioned in Section 0 we prove Lemmas IQ and E21 for the weighted version of the Hamburger 
Theorem, thereby proving Theorem 12. Ml Theorem ll.il follows as a special case of Theorem 12.81 

4.1 Proof of Lemma 12.1 L Part I 

Recall Properties 1 and 2 as well as Lemma IQ 

Property 1. The walk system contains a walk that has a self-intersection. 

Property 2. The walk system has two intersecting walks, neither of which is a 2-cycle. 

Lemma 12.11 The set of all walk systems W that satisfy either Property 1 or Property 2 can be parti¬ 
tioned into equivalence classes, each of which contributes a net zero to the permutation expansion of the 
determinant of Mh ■ 

The proof of Lemma rm is a generalization of the involution principle, the idea of which comes from 
the picture presented in Figure El Given a self-intersecting walk, changing the order of edge traversal at 
the vertex of self-intersection leads to breaking the one self-intersecting walk into two walks that intersect 
at that same vertex. Since the edge set of the collection of walks has not changed, the weight of the two 
WSP-pairs is the same. We have introduced a transposition into the sign of the permutation cycle of the 
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Figure 6: A self-intersecting cycle and its corresponding pair of intersecting cycles 


WSP-pair; this changes the sign of the WSP-pair, so these two WSP-pairs will cancel in the permutation 
expansion of the determinant of Mu- 

One can imagine that this means that to every self-intersecting WSP-pair we can associate one WSP- 
pair with two cycles intersecting. However, more than one self-intersection may occur at this same point, 
and there may be additional walks that pass through that same point. Exactly what this WSP-pair would 
cancel with is not clear. If we decide to break all the self-intersections so that we have some number N 
of walks through our vertex, it is not clear how we should sew the cycles back together. One starts to get 
the idea that we must consider all possible ways of sewing back together. Once we do just that, we have 
a family of WSP-pairs, all of the same weight, whose net contribution to the permutation expansion of 
the determinant is zero. 

This idea is conceptually simple but the proof is notationally complicated. 

If W satisfies either Property 1 or Property 2, then there is some vertex of intersection, be it either 
a self-intersection or an intersection of two walks. Our aim is to choose a well-defined first point of 
intersection at which we will build the family T. The initial term order gives an order on walks associated 
to permutation cycles; we choose the earliest walk c Xa that has some vertex of intersection. Once we 
have determined the earliest walk, we start at and follow the walk 

C Xa I Pi * Q 1 * P 2 * * * * * Pm * Qmj 

until we reach a vertex of intersection. 

In our discussion, we make the assumption that this first vertex of intersection is a vertex v* in G\. 
A similar argument exists if the first appearance occurs in Gi- Notice that at v* there may be multiple 
self-intersections or multiple intersections of walks. We will create a family T of WSP-pairs that takes 
into account each of these possibilities. 

If we want to rigorously define the breaking of a self-intersecting walk at a vertex of self-intersection, 
we need to specify many different components of the WSP-pair (C,a). First, we need to specify on 
which walk in W we are acting. Next, we need to specify the vertex of self-intersection. Since this 
self-intersection vertex may occur in multiple paths, we need to specify which two paths we interchange 
in the breaking process. 


4.2 Definitions of Breaking and Sewing 


In the following paragraphs, we define “breaking” on WSP-pairs, which takes in a WSP-pair (W,cr), one 
of cr’s permutation cycles \ a , the associated walk c Xa , paths P y and P z in c Xa , and the vertex v* in both 
P y and P z where c Xa has a self-intersection. For simplicity, we assume that v* is not a distinguished 
vertex, but the argument still holds in that case. The inverse of this operation is “sewing”. 

In this framework, c Xa has the form 

C Xc '■ Pi —* Ql P 2 ■ ■ • - 1 ” Py —*■ Qy — 1 " • • • —*■ Pz —> Qz —* • ■ • — 1 ” Pi — 1 > Ql, 


where the paths P y and P z in G\ are separated into two halves as 


Py : pfi) -> u* 


p( 2 ) 
V 
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and 


P z : P 2 (1) -> v* -> P z (2) . 

Remember that P y and P z are paths that stop over at various vertices depending on the permutation 
Xa■ The vertex v* must have adjacent stop-over vertices in each of the two paths P y and P z . Let the 
adjacent stop-over vertices in P y be v v and v v>y+1 and the adjacent stop-over vertices in P z be v lfz and 

V <Px+ 1- 

This implies Xa has the form 

Xa = (vil • ' • ‘film W 11 ’ ' ’ Vy Vy+l " ‘ Vz ‘Pz+l ' ' ' >X1 ' ' ' W\ vx ). 

We can now precisely define the result of breaking. We define xp and y 7 by splitting Xa as follows: 

Xp = Oil Py Pz+i u\ Vx ) 

and 

X 7 = (lfiy+1 ••• 

with the necessary rewriting of y 7 to have as its initial entry the value <E>(x 7 ). Define c X/3 and c Xl to be 

c x , : Pi - Qi - P2 --♦ P„ (1) Pj 2 ) - Q z -4-+ P; —»• Qi, 

and 

c X7 : Pi 1 ) - u* ^ Pf) - Q y - • ■ • -> Q z _i, 
again changing the starting vertex of c Xl to v$( Cx y 

We define the breaking of the WSP-pair with the above inputs to be the WSP-pair (W 1 , a') such that 

W' = W U {c X/3 , c Xl } \ {c Xa } 

and 

a ' = °Xa 1 XpX'y = & ■ (<Py+i Vz+i)- 

The edge set of W is equal to the edge set of W, so the weight of the modified cycle systems is the 
same as the original. Since we changed a to a' by multiplying only by a transposition, the sign of the 
modified WSP-pair is opposite to that of the original. 

4.3 Proof of Lemma 12.11 Part II 

Having defined breaking and sewing, we can continue the proof. 

For any WSP-pair (W, a) satisfying either Property 1 or Property 2, let c Xa £ W be the first walk 
in the initial term order with an vertex of intersection. Let v* be the first vertex of intersection in c y . 

A, a 

Then for all walks c with one or more self-intersections at v*, continue to break c at v* until there are no 
more self-intersections. Define the resulting WSP-pair (>V s ,cr s ) to be the simple WSP-pair associated to 
(W, a). In (W s , <t s ), there is some number N of general walks intersecting at vertex v*. There may be a 
2-cycle intersecting v* as well, but this does not matter. 

For any permutation £ £ Sjy, let £ = £ 1(2 • • • Cy b e it s cycle representation, where each Cl is a cycle. 
For each 1 < 1 < 77 , sew together walks in order: if Cl = (£ t i ■ • ■ S LS J, sew together c X5 and c xs ^ at 
v*. Sew this result together with c xs , and so on through c xs . Note that the result of these sewings is 
unique, and that every WSP-pair (W, cr) with (W s ,cr s ) as its simple WSP-pair can be obtained in this 
way, and no other WSP-pair appears. We can perform this procedure for any £ £ Sn) the sign of the 
resulting walk system is sgn(Q, a^) = sgn(£) sgn(W s , cr s ). This means that the contribution to 

the determinant of the weights of all WSP-pairs in the family T is 

T. sgn(^) sgn(W s , <J S ) wt(W s , cr s ) = sgn(W s , a s ) wt(W s , a s ) ^ sgn(^) = 0. (6) 

So elements of the same family cancel out in the determinant of Mh, giving that the contributions of all 
WSP-pairs satisfying either Property 1 or Property 2 cancel out in the hamburger determinant. • 
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4.4 Proof of Lemma 12.21 

Recall Properties 3 and 4 as well as Lemma O 


Property 3. The walk system has two intersecting walks, one of which is a 2-cycle. 

Property 4. The walk system is not minimal. 

Lemma l2.21 The set of all walk systems W that satisfy neither Property 1 nor Property 2 and that satisfy 
either Property 3 or Property 4 can be partitioned into equivalence classes, each of which contributes a 
net zero to the permutation expansion of the determinant of Mh ■ 

In the proof of Lemma E2 we do not base our family T around a singular vertex; instead, we find a 
set of violations that each member of the family has. If the WSP-pair satisfies the hypotheses of Lemma 
l2~2l then either there is some 2-cycle c : v v —► w u —► v v that intersects with some other walk or there is 
some non-minimal permutation cycle. 

Define a set of indices I C [k] of violations, of which an integer can become a member in one of two 
ways. If (W, a) is not minimal, there is at least one permutation cycle Xq with more than two consecutive 
(p’s or tc’s in its cycle notation. For any intermediary i between two (p’s or l + k between two ui’s, place 
l in I. For example, if Xa = (• • • (p' l <p" ■ ■ ■), we place i £ I. Alternatively, there may be a 2-cycle 
c : Vi —> Wk+i —*■ Vi such that either m is a vertex in some other walk c Xft or Wk+i is a vertex in some 
other cycle c X7 , or both. We declare this i to be in / as well. 

Note that any WSP-pair (W, a) satisfying either Property 3 or Property 4 has a non-empty set I. 
From our original WSP-pair, associate a minimal WSP-pair (W m , a m ) by removing any transposition \a 
from er and its corresponding 2-cycle c Xa from W, and also for any non-minimal permutation cycle xp 
we remove any intermediary (p’s or ui’s from a. We do not change the associated walk c X/3 in W since it 
still corresponds to this minimized permutation cycle. 

Let i be any element in I. Since i £ I, the 2-cycle Ci : Vi —> Wk+i —■* Vi intersects some walk of W 
either at Vi, at Wk+i , or both. So there are four cases: 


Case 1. Ci intersects a walk c Xfj 
Case 2. c, intersects a walk c x 
Case 3. Ci intersects a walk c X(3 
Case 4. Ci intersects a walk c Xfj 


at Vi and no walk at Wk+i- 
at Wk+i and no walk at Vi . 
at Vi and the same walk again at iUk+i- 
at Vi and another walk c X7 at Wk+i- 


See Figure 0 for a picture. 

In order to build the family T, we modify the minimal WSP-pair (W m , a m ) at each i with some choice 
of options, each giving a possible WSP-pair that has (W m ,cr m ) as its minimal WSP-pair. We define a 
set of two or four options Oi for each i. 

In Case 1, there are two options. Let oi be the option of including the walk Cj in W m and its 
corresponding transposition Xa in <j m . Let 02 be the option of including the intermediary <p = i in xp in 
the position where c Xfj passes through . 

[Note that we could not apply both options at the same time since the resulting y a and X/3 would not 
be disjoint and therefore is not a term in the permutation expansion of the determinant.] 

Similarly in Case 2, there are two options. Let 01 be the option of including the walk Ci in W m and 
its corresponding transposition y a in a m . Let 02 be the option of including the intermediary w = i + k 
in Xi i n the position where c x passes through Wk+i- 

In Case 3, there are four options. Let 01 be the option of including the walk c* in W m and its 
corresponding transposition Xa in er m - Let 02 be the option of including the intermediary ip = i in in 
the position where c X/3 passes through t>j. Let 03 be the option of including the intermediary ui = i + k in 
X/3 in the position where c X/3 passes through Wk+i ■ Let 04 be the option of including both intermediaries 
ip = i and u> = i + k in X /3 in the respective positions where c Xfj passes through Vi and Wk+i ■ 

In Case 4, there are four options. Let 01 be the option of including the walk q in W m and its 
corresponding transposition Xa in er m . Let 02 be the option of including the intermediary p = i in Xp in 
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1 ) 




Figure 7: The four cases in which a walk can intersect with a 2-cycle 



(1364)(25) 
sign = +1 



(12364) 
sign = -1 


Figure 8 : A family of two walk systems that cancel in the hamburger determinant 


the position where c Xl3 passes through Vi. Let 03 be the option of including the intermediary u> = i + k in 
X 7 in the position where c* passes through Wk+i ■ Let 04 be the option of including intermediary ip = i 
in X/3 in the position where c X/3 passes through V{ and intermediary ui = i + k in \ 1 in the position where 
c Xl passes through u>k+i- 

Corresponding to the associated minimal WSP-pair (W m ,<r m ) and index set /, define the family T 
to be the set of WSP-pairs (Wf, erf) by exercising all combinations of options Oj i £ Ot on (W m , cr m ) for 
each i £ I. Note that every WSP-pair derived in this fashion is a WSP-pair satisfying the hypotheses 
of Lemma l 2 ~ 2 l and is such that its associated minimal WSP-pair is (W m ,cr m ). There is also no other 
WSP-pair (W',cr') with (W m ,a m ) as its minimal WSP-pair. See Figure[E]for a canceling family of two 
walk systems, one of which is the non-minimal walk system from Figure 0 

Every WSP-pair in T has the same weight since each option changes the edge set of W by at most 
a 2 -cycle v- L — > Wk+i —> Uj, where 2 < i < k — 1 and each of those 2 -cycles contributes a multiplier of 
1 to the weight of the WSP-pair. Note that the peculiar bounds in this restriction are due to the fact 
that no permutation a can ever include two cycles through any of the vertices i’i, w^+i , Vk, or W 2 k ■ The 
sign of every WSP-pair in T is determined by the set of options applied to (W m ,cr m ). Each option Oj 
contributes a multiplicative ±1 depending on j —sgn(oi) = + 1 , sgn(c> 2 ) = — 1 , and sgn(c> 3 ) = —1 and 
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Figure 9: The symmetric difference of two matchings gives a cycle system 


sgn(o 4 ) = +1 if they exist. 

Each family T contributes a cumulative weight of 0. This is because there are the same number of 
positive (W/,£T/) € T as negative (W/, 07 ) € F. We see this since if i* £ I then any set of WSP-pairs 
that exercise the same options Oj i for all i ^ i* has either two or four members (depending on the case 
of i*), which split evenly between positive and negative WSP-pairs. Since each family contributes a net 
zero to the hamburger determinant, the total contribution from WSP-pairs satisfying the hypotheses of 
Lemma IT 2 I is zero. • 

Since every WSP-pair appears once in the permutation expansion of det Mjj, the only WSP-pairs 
that contribute to the sum are those that are simple and minimal. Therefore det Mh is the sum over 
such WSP-pairs of sgn(W, a) -wt(W, a). This establishes Theorem 12 ,.‘II and gives Theorem If. II as a special 
case. 


5 Applications of the Hamburger Theorem 

We discuss first the application of the Hamburger Theorem in the case where the region is an Aztec 
diamond, mirroring results of Brualdi and Kirkland. Then we discuss the results from the case where 
the region is an Aztec pillow, and along the way, we explain how to implement the Hamburger Theorem 
when our region is any generalized Aztec pillow. 

5.1 Domino Tilings and Digraphs 

The Hamburger Theorem applies to counting the number of domino tilings of Aztec diamonds and 
generalized Aztec pillows. To illustrate this connection, we count domino tilings of the Aztec diamond 
by counting an equivalent quantity, the number of matchings on the dual graph G of the Aztec diamond. 
We use the natural matching N of horizontal neighbors in G, as exemplified in Figure on AD 4 , as a 
reference point. Given any other matching M on G, such as in Figure |§b, its symmetric difference with 
the natural matching is a collection of cycles in the graph, such as in Figure Et- 

When we orient the edges in N from black vertices to white vertices and in M from white vertices to 
black vertices, the symmetric difference becomes a collection of directed cycles. Notice that edges in the 
upper half of G all go from left to right and the edges in the bottom half of G go from right to left. 

Since the edges of N always appear in the cycles, we can contract the edges of N to points; then the 
graph retains its structure in terms of the cycle systems it produces. This new graph H is of the form 
in Figure UHL This argument shows that the number of domino tilings of an Aztec diamond equals the 
number of cycle systems of this new condensed graph, called the region’s digraph. 

We wish to concretize this notion of a digraph of the Aztec diamond AD n . Given the natural tiling 
of an Aztec diamond consisting solely of horizontal dominoes, we place a vertex in the center of every 
domino. The edges of this digraph fall into three families. From every vertex in the top half of the 
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Figure 10: The hamburger graph for an Aztec diamond and an Aztec pillow 


diamond, create edges to the east, to the northeast, and to the southeast whenever there is a vertex 
there. From every vertex in the bottom half of the diamond, form edges to the west, to the southwest, 
and to the northwest whenever there is a vertex there. Additionally, label the bottom vertices in the top 
half v% through v n from west to east and the top vertices in the bottom half w n +i through W2n > also west 
to east. For all i between 1 and n, create a directed edge from u ?: to w n+ i and from w n+ i to Vj. The 
result when this construction is applied to AD 4 is a graph of the form in Figure ITOl i. By construction, 
this graph is a hamburger graph. 

Theorem 5.1. The digraph of an Aztec diamond is a hamburger graph. 

Since both the upper half of the digraph and the lower half of the digraph are strongly planar, there 
are no negative cycle systems. This implies that the determinant of the corresponding hamburger matrix 
counts exactly the number of cycle systems in the digraph. 

Corollary 5.2. The number of domino tilings of an Aztec diamond is the determinant of its hamburger 
matrix. 

5.2 Applying the Hamburger Theorem to Aztec Diamonds 

In order to apply Theorem I I .1 1 to count the number of tilings of AD n , we need to find the number of 
paths in the upper half of D from v t to Vj and the number of paths in the lower half of D from w n +j 
to w n +i. The key observation is that by the equivalence in Figure eh we are in effect counting the 
number of paths from (i,i) to (j. j) using steps of size (0,1), (1,0), or (1,1) that do not pass above the 
line y = x. This is a combinatorial interpretation for the (j — *)-th large Schroder number. The large 
Schroder numbers So,..., S 5 are 1,2,6, 22,90, 394. This is sequence number A006318 in the Encyclopedia 
of Integer Sequences EH 

Corollary 5.3. The number of domino tilings of the Aztec diamond AD n is equal to 

#AD n = det 

where S n is an upper triangular matrix with the i-th 


_ r qT > 

1 n <-'n 

large Schroder number on its ith superdiagonal. 
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lattice paths in the first quadrant 


394 
90 
22 
6 
2 
1 

Brualdi and Kirkland prove a similar determinant formula for the number of tilings of an Aztec 
diamond in a matrix-theoretical fashion based on the n{n + 1) x n(n + 1) Kasteleyn matrix of the graph 
H and a Schur complement calculation. The Hamburger Theorem gives a purely combinatorial way to 
reduce the calculation of the number of tilings of an Aztec diamond to the calculation of a 2 n x 2 n 
Hamburger determinant. 

Following cues from Brualdi and Kirkland, we can reduce this to an n x n determinant via a Schur 
complement. For uses and the history of the Schur Complement, see the new book 1X3 by Zhang. In the 
case of the block matrix Mjj in Equation £3, taking the Schur complement of B in Mu gives 

det Mh = det B ■ det (A + DiB~ 1 D 2 ) = det (A + DiB~ 1 D 2 ), (7) 


Figure 11: The equivalence between paths in D and 


For example, when n = 6, the matrix Sq is 


Se = 


1 2 6 22 90 

0 1 2 6 22 

0 0 12 6 
0 0 0 1 2 

0 0 0 0 1 

0 0 0 0 0 


since B is a lower triangular matrix with l’s on the diagonal. In this way, every hamburger determinant 
can be reduced to a smaller determinant of a Schur complement matrix. In the case of a simple hamburger 
graph where D 2 = D\ = /, the determinant reduces further to det(A + 5 _1 ). Lastly, in the case where 
the hamburger graph is rotationally symmetric, B = JAJ where J is the exchange matrix, which consists 
of l’s down the skew-diagonal and 0’s elsewhere. This implies that we can write the determinant in terms 
of just the submatrix A , i.e., det (A + JA _1 J). (Note that J -1 = J.) 


Corollary 5.4. The number of domino tilings of the Aztec diamond AD n is equal to det (S n + J n S n 1 J n ), 
where J n is the n x n exchange matrix. 


In the case of a hamburger graph H , we call this Schur complement, A + S , a reduced hamburger 
matrix. In the Aztec diamond graph example above, we can thus calculate the number of tilings of the 
Aztec diamond AD 6 as follows. The inverse of Se is 



1 

-2 

—2 

—6 

-22 

-90 

0 

1 

-2 

-2 

-6 

-22 

0 

0 

1 

-2 

-2 

-6 

0 

0 

0 

1 

-2 

-2 

0 

0 

0 

0 

1 

-2 

0 

0 

0 

0 

0 

1 
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Figure 12: The digraph of a generalized Aztec pillow from the digraph of an Aztec diamond 


which implies that the determinant of the reduced hamburger matrix 


2 

2 

6 

22 

90 

394 

-2 

2 

2 

6 

22 

90 

-2 

-2 

2 

2 

6 

22 

-6 

-2 

-2 

2 

2 

6 

-22 

-6 

-2 

-2 

2 

2 

-90 

-22 

-6 

-2 

-2 

2 


gives the number of tilings of ADq . 

Brualdi and Kirkland were the first to find such a determinantal formula for the number of tilings of 
an Aztec Diamond [2j. Their matrix was different only in the fact that each entry was multiplied by (—1) 
and there was a multiplicative factor of (—l) n . Brualdi and Kirkland were able to calculate the sequence 
of determinants, {det M„}, using a J-fraction expansion, which only works when matrices are Toeplitz 
or Hankel. 

5.3 Applying the Hamburger Theorem to Aztec Pillows 

A generalized Aztec pillow can be produced by restricting the placement of certain dominoes in an Aztec 
diamond in order to generate the desired boundary. We define the digraph of a generalized Aztec pillow 
to be the restriction of the digraph of an Aztec diamond to the vertices on the interior of the pillow. An 
example is presented Figure El Since the generalized Aztec pillow’s digraph is a restriction of the Aztec 
diamond’s digraph, we have the following corollary. 

Corollary 5.5. The digraph of an Aztec pillow or a generalized Aztec pillow is a hamburger graph. 

Aztec pillows were introduced in part because of an intriguing conjecture for the number of their 
tilings given by Propp [El- 

Conjecture 5.6 (Propp’s Conjecture). The number of tilings of an Aztec pillow AP n is a larger 
number squared times a smaller number. We write ffAP n = l^s n . In addition, depending on the parity 
ofn, the smaller number s n satisfies a simple generating function. For AP^m, the generating function is 

OO 

S2mX m = (1 + 3x + x 2 — £ 3 )/(l — 2x — 2x 2 — 2x 3 + X 4 ). 

m =0 

For AP 2 m+i, the generating function is 

OO 

S2m+lX m = (2 + x + 2x 2 — x 3 )/{l — 2x — 2x 2 — 2x 3 + x 4 ). 

m —0 
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Figure 13: The equivalence between paths in D and lattice paths in the first quadrant 


My hope for the Hamburger Theorem was that it would give a proof of Propp’s Conjecture. We shall see 
that although we achieve a faster determinantal method to calculate #AP ra , the sequence of determinants 
is not calculable by known methods. 

Using the same method as for Aztec diamonds, creating the hamburger graph H for an Aztec pillow 
gives Figure Eb- Counting the number of paths from Vi to Vj and from w n +j to w n +i in successively 
larger Aztec pillows gives us the infinite upper-triangular array S = ( Sij ) of modified Schroder numbers 
defined by the following combinatorial interpretation. Let Sij be the number of paths from ( i,i ) to (j,j) 
using steps of size (0,1), (1, 0), and (1,1), not passing above the line y = x nor below the line y = x/ 2 . 
The equivalence between these paths and paths in the upper half of the hamburger graph is shown in 
Figure El The principal 7x7 minor matrix S 7 of S is 
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0 

0 
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0 

0 

0 

0 

0 
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0 

0 

0 

0 
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This table of numbers is sequence number Al 14292 in the Encyclopedia of Integer Sequences El 
Theorem 5.7. The number of domino tilings of an Aztec pillow of order n is equal to 


#AP n = det 


S n In 

In Jn Sn Jn 


where S n is the n x n principal submatrix of S and J n is the n x n exchange matrix. 

As in the case of Aztec diamonds, we can calculate the reduced hamburger matrix through a Schur 
calculation. The inverse of S 6 is 


S, 


-l 


and the resulting reduced hamburger 


M e = 



' 1 

-1 

0 

0 

0 

0 


0 

1 

-2 

-1 

-2 

-5 


0 

0 

1 

-2 

-2 

-5 


0 

0 

0 

1 

-2 

-2 


0 

0 

0 

0 

1 

-2 


0 

0 

0 

0 

0 

1 

matrix 

for AP 6 is 




- 

2 

1 

2 

5 

16 

57 

-2 

2 

2 

5 

16 

57 

-2 

-2 

2 

2 

6 

21 

-5 

-2 

-2 

2 

2 

6 

-5 

-2 

-1 

-2 

2 

2 

0 

0 

0 

0 

-1 

2 


This gives us a much faster way to calculate the number of domino tilings of an Aztec pillow than was 
known previously. We have reduced the calculation of the 0(n 2 ) x 0(n 2 ) Kasteleyn-Percus determinant 
to an n x n reduced hamburger matrix. To be fair, the Kasteleyn-Percus matrix has —1,0, and +1 entries 
while the reduced hamburger matrix may have very large entries, which makes theoretical running time 
comparisons difficult. Experimentally, when calculating the number of domino tilings of AP 14 using 
Maple 8.0 on a 447 MHz Pentium III processor, the determinant of the 112 x 112 Kasteleyn-Percus 
matrix took 25.3 seconds while the determinant of the 14 x 14 reduced hamburger matrix took less than 
0.1 seconds. 

Whereas we have a very understandable determinantal formula for the number of tilings of the region, 
this does not translate into a proof of Propp’s Conjecture because we cannot calculate the determinant 
of the matrices M n explicitly. We cannot apply a J-fraction expansion as Brualdi and Kirkland did since 
the reduced hamburger matrix is not Toeplitz or Hankel. 


6 Generalizations 

In this section we present a counterexample to a natural generalization of the Hamburger Theorem and 
discuss an intriguing generalization of Propp’s Conjecture. 
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Figure 14: A simple generalized hamburger graph 
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Figure 15: The ten cycle systems for the generalized hamburger graph in Figure ITU 


6.1 Counterexample to Generalizations of the Hamburger Theorem 

The structure of the hamburger graph presented in Figure ^ seems restrictive, so the question naturally 
arises whether it is somehow necessary. Can the edge set E 3 between graphs Gi and Gi be an arbitrary 
bipartite graph? The answer in general is no. 

Take for example the simple graph H in Figure As one can count by hand, there are 10 distinct 
cycle systems in H, as enumerated in Figure ITT)! Creating the hamburger matrix that would correspond 
to this graph gives 

'1 1 110 0 " 

0 1 10 10 

0 0 1 0 0 1 

Mh ~ -1 0 0 1 0 0 ’ 

0 0-1110 
0-10111 

where the lower left block describes paths from uq, uq, and uq to iq, iq, and tq and we have changed 
no other block of the matrix. This seems the most logical extension of the hamburger matrix. 

The determinant of this matrix is —5. Even though one might expect to need a new sign convention on 
cycle systems in generalized hamburger graphs, any sign convention would necessarily conserve the parity 
of the number of cycle systems. Therefore there is no +1/—1 labeling of the cycle systems in FigurelTHlthat 
would allow clet Mu = X^ c eC s § n ( c )’ This means that either the matrix for this generalized hamburger 
graph is not correct or that the theorem does not hold in general. 

The only hope for my idea of how to generalize the hamburger matrix is what happens when one 
considers the difference between even and odd permutations. In the counterexample above, the edge set 
E 3 sends vertex uq to vertex tq, uq to U3, and uq to tq. This can be thought of as an odd permutation 
of a hamburger’s edge set (sending uq to tq, uq to tq, and uq to V3). In my limited calculations so far, 
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q = 3 


q = 5 




q = 7 


q = 9 


Figure 16: Plots of values of s ntq /s n - 2 , q for q = 3, 5, 7,9 


whenever the permutation of the natural edge set is even, the parity of the number of cycle systems, 
counted by hand, is equal to the parity of the determinant of the matrix. I have not found a sign 
convention for these graphs that explains the calculated determinant, but being equal in parity allows 
for hope of a possible extension to the Hamburger Theorem. 

6.2 Generalizing Propp’s Conjecture 

A key motivation in the study of Aztec pillows is Propp’s Conjecture for 3-pillows, presented in Section 
PI Through experimental calculations, I have expanded and extended the conjecture. 

For one, I extend the conjecture to odd pillows. For any odd number g, we can use the notation AP% 
to represent the n-th g-pillow for odd q. Recall from the introduction that this is a centrally symmetric 
region with steps of length q and central belt of size 2 x 2 n. Calculating the number of tilings of AP%, 
AP%, and AP% for n at least up to 40 gives strong evidence that #AP® for q odd is always a larger 
number squared times a smaller number. As with 3-pillows, we can write #AP% = 1^ q s nA , where s„^ q is 
relatively small but is not necessarily the squarefree part of #AP %. 

Additional structure in these smaller numbers also appears. Propp’s Conjecture gives an explicit 
recurrence for the smaller-number values when q = 3. In the case q = 5, there is no linear, constant 
coefficient fc-th degree recurrence for any k up to 20. However in each case (g = 5, g = 7, and g = 9) the 
structure is undeniable. Plotting the values of s n , q /s n - 2 ,q gives a remarkable damped sinusoidal graph 
converging to some number p. (See Figure UTil) 

In addition, plotting the values of s n , q /s n - i, g depending on whether n is even or odd gives two damped 
sinusoidal graphs converging to some values r and 2r respectively. This implies p = 2r 2 . Approximate 
values for p are given in Tabled We have included the values for Aztec diamonds (g = 1) and the original 
Aztec pillows (q = 3). 

Even more appears to be true. By construction, the dual graphs to all odd pillows are 2-even- 
symmetric—that is, centrally symmetric so that antipodes are an even path length apart. Jockusch p] 
proved that the number of matchings of such graphs is the sum of two squares. 
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q Approximate limit of s n ^ q /s n - 2 , q 
1 2 “ 
3 2.890053636 

5 3.0821372 

7 3.145 

9 3.18 


Table 1: p values for various values of q 


Calculating these squares gives us a better understanding of the larger numbers l n appearing in 
Propp’s Conjecture, and the larger numbers l ntq in general. In the case of odd pillows, we can write 
#AP% = q for explicitly known positive integers a n , q and b n ^ q . Surprisingly, in each case I 

considered, l nA divides both a n , q and This implies that the smaller number is also a sum of 

squares, which may lend some insight into its recurrence formula. 

This allows us to formulate a new and improved conjecture about odd pillows. 

Conjecture 6.1. The number of tilings of AP% is a larger number squared times a smaller number. 
There exist values and such that we can write ffAP% = l^^n^. The value s n satisfies the 
following structure. The ratio of S 2 n+\, q / S 2 n,<? to S 2 n+ 2 ,q/S 2 n+i,q is exactly 2 in the limit. In addition, 
for the explicit values a nA and b nA from ffAPif’s sum of squares representation, the value l n ^ q from above 
divides both a n , q and b n ^ q . 

6.3 Conclusion 

The Hamburger Theorem allows us to calculate the number of cycle systems in generalized Aztec Pillow 
graphs more quickly. It also allows us to count the number of cycle systems in many previously inaccessible 
graphs, such as non-planar hamburger graphs. Just as Gessel and Viennot’s result has found many 
applications, I hope the Hamburger Theorem will be useful as well. 

6.4 Remarks 

This article is based on work in the author’s doctoral dissertation. A preliminary version of this arti¬ 
cle appeared as a poster at the 2005 International Conference on Formal Power Series and Algebraic 
Combinatorics in Taormina, Italy. The author would like to thank Henry Cohn and Tom Zaslavsky for 
numerous intriguing discussions and useful corrections. 
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